Infinite-randomness quantum Ising critical fixed points 
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We examine tlie ground state of the random quantum Ising model in a transverse field using a gen- 
eralization of the Ma-Dasgupta-Hu renormalization group (RG) scheme. For spatial dimensionality 
d — 2, we find that at strong randomness the RG flow for the quantum critical point is towards 
an infinite-randomness fixed point, as in one-dimension. This is consistent with the results of a 
recent quantum Monte Carlo study by Pich, et al., including estimates of the critical exponents 
from our RG that agree well with those from the quantum Monte Carlo. The same qualitative 
behavior appears to occur for three-dimensions; we have not yet been able to determine whether or 
not it persists to arbitrarily high d. Some consequences of the infinite-randomness fixed point for 
the quantum critical scaling behavior are discussed. Because frustration is irrelevant in the infinite- 
randomness limit, the same fixed point should govern both ferromagnetic and spin-glass quantum 
critical points. This RG maps the random quantum Ising model with strong disorder onto a novel 
type of percolation/aggregation process. 



I. INTRODUCTION 

Systems with quenched randomness and many degrees 
of freedom may be divided into three classes based on 
their coarse-grained behavior in the low-energy, low- 
frequency and/or long-distance limit. First, there are 
many systems where the quenched disorder is irrele- 
vant in the renormalization-group sense. Such systems, 
even though they are spatially inhomogcncous at the 
microscopic scale, become asymptotically homogeneous 
at macroscopic scales; their coarse-grained, low-energy 
behavior is the same as some "pure" system without 
quenched disorder. In this case, the renormalization- 
group fixed point governing the coarse-grained system 
is at zero quenched randomness. The second possibil- 
ity is systems controlled by fixed points with nonzero, 
but finite, quenched randomness. In this case the coarse- 
grained behavior is spatially inhomogeneous, but the rel- 
ative magnitude of the inhomogeneities remains finite at 
the fixed point. Examples of this second class include 
spin glasses and other glassy phases, as well as vari- 
ous critical points with randomness. The third possibil- 
ity, which is the subject of this paper, occurs when the 
quenched randomness and thus the relative magnitude 
of the inhomogeneities grows without limit as the sys- 
tem is coarse-grained. So far, we know of a few infinite- 
randomness fixed points which comprise this third class; 
all of them are one-dimensional (d = 1) quantum ground 
states. These are the random singlet i^fptes of certain 
random antiferromagnetic spin chains ,Qa the quantum 
critical uoint of the random quantum Ising (and Potts) 
chain,QLl and quantum critical points separating raru 
dom singlet states and the Ising antiferromagnetic phasetl 
or the-.ljlaldane state in the random spin-1 Heisenberg 
chain .Era 

A natural question is whether such infinite-randomness 
fixed points can govern the behavior of physical systems 



with spatial dimensionality d > 2. 

Here we study the simplest random quantum system 
that can exhibit a phase transition, the random ferro- 
magnetic Ising model in a transverse field, which, in all 
dimensions, has a quantum critical point at zero temper- 
ature. We focus_pn the strong-randomnesa^regime, using 
a generalizationfl of the Ma-Dasgupta-HuO renormaliza- 
tion group (RG) scheme. This approximate RG is exact 
in the limit of infinite randomness and thus can in princi- 
ple yield exact results for the scaling behaviot,of systems 
governed by infinite randomness fixed points.B In one di- 
mension, this RG has been used to analyze the various 
infinite-randomness fixed points mentioned above and 
many results can be obtained analytically, in particular 
for the random quantum Ising model. In higher dimen- 
sions, the renormalization group cannot (to our knowl- 
edge) be carried out analytically; in this paper we analyze 
its general structure and consequences, implement it nu- 
merically and examine a simple approximation to it. For 
d = 2 and 3, we find that the renormalization group flow 
on the critical manifold for strong randomness is indeed 
towards even stronger randomness, as in one-dimension, 
indicating that the quantum critical behavior is governed 
by the infinite-randomness critical fixed point. For d — 2 
we have studied the RG flow thoroughly enough to be 
fully confident that this is the case, and it appears to 
remain true for d — 3. ■— ■ 

For d < 4, the Harris criterionliJ indicates that at the 
pure Ising quantum critical point weak randomness is rel- 
evant, with the RG flow towards stronger randomness; 
thus the simplest scenario for two- and three-dimensions 
is that the inflnite-randomness fixed point governs the 
random quantum critical point with any randomness, as 
in one-dimension. A recent quantum Monte Carlo study 
by Pich, et al. of the ferromagnetic model with moder- 
ate randopmess in two-dimensions is consistent with this 
pictureO 

What do we mean by infinite randomness? This 



meanga'B that as the system is coarse-grained and the 
characteristic energy scale decreases, the distributions 
of the logarithms of the magnitudes of the terms in the 
renormalizcd hamiltonian become arbitrarily broad. As 
a result, the ratio of the magnitude of any two terms 
approaches either zero or infinity. In this limit each 
renormalized coupling completely dominates any weaker 
renormalizcd couplins, so even the Ising spin glass be- 
comes unfrustratedJlj Thus we expect that the same 
infinite-randomness fixed point governs both the random 
ferromagnetic and the spin-glass quantum critical points, 
with simple modifications to account for the antifcrro- 
magnetic bonds in the latter case. We should note that 
recent quantum Monte Carlo studies of the spin-glass 
case had concluded that for two and three dimensions tba 
scaling near the quantum critical point is conventional,£3 
implying a finite-randomness fixed point, in contrast to 
what we are proposing here. But these finite-size scaling 
Monte Carlo studies looked at a rather small size range 
and did not look at distributions of physical properties; 
thus they were not sensitive to the scaling towards infi- 
nite randomness that we now believe occurs for d = 2 
and probably for higher d. 



II. CLUSTER RG 

We will study the quantum spin-I/2 Ising model with 
random ferromagnetic interactions, positive transverse 
fields hi, and moments fif. 



Jjk ~max(Jjfc,— -— ). 
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where the ct's are the Pauli spin matrices. The interac- 
tions Jij dictate the geometry and thus the dimensional- 
ity of the system, with the strong interactions being pre- 
dominantly between nearby pairs of spins. The uniform 
ordering field H — H^ is zero or small; it is included only 
to probe the system's magnetization and susceptibilities. 
The cluster RG finds the system's ground state by suc- 
cessively eliminating the highest energy degrees of free- 
dom. At each step, we find the largest term in the hamil- 
tonian, which is either a transverse field or an interac- 
tion; its strength 17 sets the (maximum remaining) energy 
scale. If the largest term is the field on spin i, SI = hi, 
that spin is put in the ground state, erf — I, of the lo- 
cal field term and virtual excitations to the other state 
(af = — 1) are treated in second-order perturbation the- 
ory. For this to be a valid approximation, the field must 
be much stronger than all this spin's interactions, which 
is true only in the strong-randomness limit. This step 
eliminates the spin i, and generates new interactions of 
the form J'^ « JjiJik/hi. But some of these new inter- 
actions may be negligible compared to the interactions 
that were already present, so the full renormalized inter- 
actions are given, for each pair (j, k), by 



We use the maximum here because in the strong- 
randomness limit the sum of two nonnegative numbers 
of very different magnitudes is well approximated by the 
larger number. 

If, on the other hand, the largest term in the hamilto- 
nian is an interaction, fJ — Jij, the two spins involved are 
combined into one new spin — a cluster — whose two states 
represent the two ground states {al = o-j = ±1) of that 
interaction. Again, the virtual excitations to the states 
that are eliminated are treated in second-order perturba- 
tion theory. With the new spin labeled i (an arbitrary 
choice) the renormalized field is 



h' 



hihj 
'J a 



and the renormalized moment is simply 

m'j =/At+ Mi- 



es) 



(4) 



The renormalized interactions are, for each remaining fc. 



J-fc « max(Jife, Jjfe). 



(5) 



The net result in both cases is the elimination of one 
spin along with the various renormalizations and recon- 
nections of the lattice. The major complication for d > \ 
is that the RG does not preserve the lattice structure 
but instead generates a disordered and strongly corre- 
lated network of sites — i.e. clusters — and bonds connect- 
ing them. 

The action of the cluster RG is a novel aggrega- 
tion/annihilation process. When the strongest term is a 
field, the corresponding cluster is removed (annihilated), 
while when it is an interaction the two clusters that it 
connects are aggregated into one cluster. The clusters 
thus represent sets of the original spins that are strongly 
correlated. In the paramagnetic phase, the annihilation 
process dominates and no large clusters are formed, while 
in the ordered phase the aggregation dominates produc- 
ing arbitrarily large clusters. For d > 1, in the ordered 
phase an infinite percolating cluster appears at a finite 
energy scale. When this occurs, the infinite cluster is 
represented by a single renormalized spin that has an in- 
finite coordination number and an infinite magnetic mo- 
ment. Thus the topology of the network has by this en- 
ergy scale completely changed from that of the "bare" 
d-dimensional lattice which had only short-range inter- 
actions. This emphasizes that this RG is not simply a 
"real-space" RG; it is more precisely an "energy-space" 
RG. 

The zero-temperature quantum critical point is a new 
type of percolation transition at which an infinite cluster 
first appears (at the quantum critical point it only ap- 
pears in the. zero energy limit) , as pointed out by Mon- 
thus, et aln. The quantum critical point occurs when 



the annihilation and aggregation processes balance, so 
that arbitrarily large clusters are produced but no single 
cluster dominates at any finite energy scale. 

Because of the multiplicative structure of the RG re- 
cursion relations, it is convenient to write the renormal- 
ized interactions after decimating down to energy scale 
fi in logarithmic variables as: 



hi — U,e 



-ft 



(6) 

(7) 



The RG equations involve the full joint distribution of 
the positive quantities {A} and {Cy}- The single-field 
distribution function R{(3) and the single-bond distribu- 
tion function P(C) provide partial but important infor- 
mation about the full joint distribution of all the cou- 
plings. In general, the coordination number is not fixed, 
and if there are N spins (clusters) remaining, there may 
be up to N{N — l)/2 interactions. However, only of or- 
der TV of the strongest interactions matter because the 
weaker ones will be overruled by the stronger ones in the 
action of the RG for strong randomness. Thus we nor- 
malize the bond distribution function P(C) per remaining 
spin, so that the total number of bonds in the system is 
TV/q P{C)dC, and for d > 1 we define the "width", w, 
of this bond distribution by 



/ PiOdC^i. 

Jo 



(8) 



This width thus includes only the N strongest bonds. 
[Note that in one dimension, there are always exactly N 
nearest neighbor bonds and this normalization of P(C) 
coincides with that of Ref.0. If further neighbor bonds ex- 
ist in one dimension, they will be irrelevant at low energy 
scales and can be ignored.] The condition for validity of 
our approximate RG is that the widths of the distribu- 
tions -R(/3) and P(C) be large; for the RG to be asymp- 
totically exact, these widths should tend to infinity as 
r2 ^ 0; this is the indication that the RG flow goes to an 
infinite-randomness fixed point. We now discuss the gen- 
eral structure of such a putative infinitc-randomness crit- 
ical fixed point. We will assume that the simplest scaling 
scenario for an_ijifinite-randomness fixed point occurs in 
our cluster RGtJ — this is the case in one-dimension,Q and 
is consistent with our numerics for two-dimensions. 



III. CRITICAL FIXED POINT 

We first consider the behavior at the quantum criti- 
cal point. We expect that at the critical fixed point the 
distributions i?(/3; F) and P(C; T) will asymptotically be 
given by the simple scaling forms 



i?(/3; T)df3 = B{l3/T)dl3/T , 

PiC;r)dC^z{c/T)dC/r 



(9) 
(10) 



for large F; where 



In 



n 



>o 



(11) 



is the logarithm of the energy scale, relative to an energy 
scale flo set by the properties of the bare hamiltonian. 
As F is increased by dF, the fractional decrease in the 
number of spins or clusters is {Zq + Bo)dr/r with 



Bo = 5(0), Zo = Z{0). 



(12) 



The density of remaining clusters per unit volume thus 
decreases under renormalization as 



n-r 



(13) 



This gives the basic relationship between the length scale, 
L, and the energy scale at the quantum critical point: 



F = ln 



n 



i" 



(14) 



with ijj — d/{Zo + Bq) < 1. Note that this is very dif- 
ferent from conventional power-law scaling; here it is the 
logarithm of the energy scale that varies as a power of 
the length scale. Since this is associated with the "tun- 
neling" events by which clustpjs flip, it has been dubbed 
"tunneling dynamic scaling" .Lj 

The typical moment of a cluster — the number of 
strongly correlated "active" spins in it (i.e., those not 
yet decimated) — scales as 



M 



p0 



(15) 



so that the fractal dimension of the set of active spins in 
a cluster is 



df = (flip. 



(16) 



This determines the decay of the average spin-spin cor- 
relation function at the critical point. 



A. Critical Correlations 



The correlation function between two spins at distance 



r, 



Gij =< erf o-| > 



(17) 



is a random quantity with a very broad probability dis- 
tribution for large r = jr^ — rj|. We first consider the 
correlation function of a typical pair of spins with sepa- 
ration r; typical spin pairs are never active in the same 
cluster and have only weak correlations that fall off, at 
criticality, as a stretched exponential function of distance. 



lnGtyp(r)~r'>, 



(18) 



with tfj < I. These correlations arise from the lowest 
order perturbative corrections to the decimation of spin 



clustersErcil. When a cluster is decimated at energy scale 
51, each of the effective spins on its neighboring clusters — 
in more conventional terms the perturbatively modified 
wavefunctions that are labeled by the remaining effective 
spins — will acquire a component of the decimated clus- 
ter's spin whose magnitude is of order J/fi, with J the 
effective coupling that links the neighboring cluster to 
the decimated cluster. Likewise, when these neighboring 
clusters are decimated, an even smaller component of the 
original spin will be acquired by the remaining clusters. 
Correlations between two spins i and j which are never 
active in the same cluster thus occur when a surviving 
cluster contains simultaneously a component of both the 
spins i and j. The correlation function Gij is then deter- 
mined by the maximum over all such mutual clusters that 
occur, at any energy scale, of the product of the the two 
spins' components contained by the mutual cluster. Typ- 
ically, the smallest of the multiplied perturbative factors 
that determine In G will dominate; these arise from the 
stage at which the two spins first have a component on 
a mutual cluster. Since this occurs when the remaining 
cluster sizes are of order the separation r between the 
spins of interest, the dominant perturbative factor will 
be of order e~^'J^ with F ~ r'^ and Kjj jandom and of 
order one. This yields the result Eq.@.&y 

We should note that there is another mechanism by 
which spins become correlated. At any step of the renor- 
malization, when a cluster with field h is decimated, 
higher order perturbative effects will give components of 
the decimated spins on all the remaining clusters, with 
magnitudes that involve products over all the bonds con- 
necting the decimated cluster to the remaining clusters, 
of factors of the form Jjk/h. These by themselves would 
give rise to simple exponential decay of typical correla- 
tions as occurs in conventional disordered phases. In con- 
trast to these, the contributions to the correlations that 
come from the cumulative effects of the successive lowest 
order perturbative terms discussed above will have simi- 
lar form but with each of the {Jjk} being divided by an 
effective field from one of the later stages of the RG; these 
are always smaller than h. Thus the cumulative lowest 
order contributions discussed above will always dominate 
over the simple exponential decays from the higher order 
perturbative effects. We hence conclude that the ip from 
the RG must be less than one. 

The average correlation function G{r) behaves quite 
differently than typical correlations. It is dominated by 
the rare spin pairs which are active in the same cluster at 
some energy scale; such pairs of spins have correlations of 
order one, reduced from one only by the short scale high 
energy fluctuations which arc not included in the approx- 
imate RG. As a result, the average correlation function 
is proportional to the probability of the two spins being 
active in the same cluster at some energy scale. At crit- 
icality this occurs — if at all — at scale F ~ r'f' , and the 
resulting average correlation function hence falls off as a 
power law: 



G{r) 



-2(d-0V) 



(19) 



This is an example of the radically different scaling be- 
havior of the typical and average quantities _tjhat is one 
hallmark of infinite-randomness fixed points. tHij 

Thermodynamic properties involve averaging over the 
whole system and will hence be dominated, as are the av- 
erage correlations, by rare clusters. The low temperature 
susceptibility to a small ordering field H {H ^ T) can be 
found easily by stopping the RG at energy scale Q ^ T. 
For small T, almost all the decimated spins are frozen 
and hence non-magnetic, while almost all the remaining 
clusters have effective transverse fields and interactions 
between them which are much less than T. They are 
hence essentially free and have independent Curie sus- 
ceptibilities yielding an overall susceptibility at low tem- 
peratures near the quantum critical point of 



X 



llnTp'^"'^/'^ 



T 



(20) 



The magnetization in a small ordering field H (at 
T <^ H) can be found similarly: the RG is stopped when 
the typical magnetic energy, iJ/i, of a cluster is of or- 
der O. The decimated spins are non-magnetic while the 
remaining clusters are almost perfectly polarized by the 
field. This yields, at the quantum critical point, a mag- 
netization proportional to the fraction of spins which are 
still active. 
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B. Simple Approximate RG 

In order to get a better feeling for the scaling be- 
havior, it is useful to study a simple approximation to 
the RG flows which is exact in one-dimension and in 
some respects good for higher d; this consists of ig- 
noring correlations among the fields and between the 
fields and the bonds, but allowing correlations among 
the bonds. In this approximation, the evolution equa- 
tion for the field distribution R{f3; F) depends on the 
bond distribution only through P(0; Fl||-and is identical 
to its one-dimensional form (see Refs.BEJ). Solving for 
the fixed point gives scaling distributions with Zq = 1 
and B{9) — Boe~^°^ with Bq undetermined (the exact 
solution for one-dimension has Bq — 1). Within this 
approximation the exponents are tp — d/ {I + Bq) and 
(j) = (1 + VI + 4i3o)/2. Our numerical studies of the 
RG flows in two-dimensions show that the log-field dis- 
tribution is very close to the simple exponential form, but 
the estimated critical exponents for two-dimensions differ 
somewhat from this simple approximation; this must be 
due to correlations among the fields and between them 
and the bonds. 



IV. NUMERICAL RG STUDY OF THE CRITICAL 
FIXED POINT IN TWO-DIMENSIONS 

In order to do better and certainly to test the conjec- 
ture of a controlling infinite-randomness fixed point, we 
must implement the strong-randomness cluster RG nu- 
merically. The formulation is the same for a general ran- 
dom network, but we are of course interested in systems 
that can arise from finite-dimensional lattices with short 
range interactions. We have thus studied the RG flows 
with "initial conditions" of finite d-dimensional lattices 
for two and three dimensions. The program has been 
tested by verifying that it reproduces (within statistical 
errors) the analytical results for one-dimension. 

For rf > 1, many weak interactions are generated that, 
for the larger lattices studied, can not all be stored. Be- 
cause of this, we keep only interactions above a mini- 
mum strength Jmin, the smallness of Jmin being limited 
by computer memory capacity and speed. Since the dis- 
carded bonds (those with J < Jmin) could not have gener- 
ated stronger bonds, for the renormalization down to any 
energy scale with D, > Jmin the RG decimation sequence 
is not affected at all, and all the bonds with Jmin < J <^ 
are retained; thus we know all the fields and all the bonds 
with < C < Cm = ln(r^/Jmin)- However, under renor- 
malization Cm decreases and F = ln(rio/^) increases so 
that the range of C/F in the scaled distribution Z(C/F) 
that we can study steadily decreases as the system is 
coarse-grained. We have been quite conservative and do 
not look at all at the "contaminated" low-energy part of 
the bond distribution (J < Jmin)- The limits on memory 
are most restrictive at the earliest stages of the RG dec- 
imation, where the number of clusters is large, so this is 
when Jmin must be set the largest. It might be possible 
to let Jmin decrease later in the decimation and recover 
more of the renormalized bond distribution with control- 
lable errors. We have not explored this possibility. 

We start with initial conditions of systems of up to 10^ 
spins with random short-range interactions and random 
transverse fields independently chosen from specified ini- 
tial probability distributions. We run up to 1000 samples 
for each initial probability distribution to reduce the sta- 
tistical errors. For each sample we measure properties of 
the system when the energy scale passes (under renormal- 
ization) a predefined set of energiesEd, and then average 
these properties over different samples. 

To reduce transients as much as possible, the shapes 
of the initial distributions R{/3) and P(C) are chosen to 
approximate, as best as we can, the renormalized criti- 
cal point distributions which we observe. However, these 
initial conditions are missing any correlations among the 
fields and interactions which certainly exist in the full 
joint distribution at the critical fixed point. Thus when 
we run the RG it does show a fairly strong transient be- 
havior as these correlations are generated and the fixed 
point is approached. So far, we have only a limited un- 
derstanding of these transients and the correlations that 



are generated and we do not have a systematic way of 
controlling them; we do, however, monitor the simplest 
types of correlations and they do appear to stabilize after 
the initial transient in the RG. 

In our numerics we primarily concentrate on the indi- 
vidual field and bond distributions i?(/3; F) and P(C;F); 
these are partial but significant indicators of what is hap- 
pening in the system's full joint probability distribution. 
For two-dimensions we find that at the critical point both 
distributions do become broader under the action of the 
RG and the flow towards stronger randomness is clear. 
This flow is weaker but nevertheless is clearly apparent 
for three-dimensions also. 

Numerically we find that under the RG the field distri- 
bution maintains fairly accurately a simple exponential 
form, 



-Ro{r)p 



with 



i?(/3;F)^i?o(r)e 



i?o(r) = i?(/3 = 0;F). 



(22) 



(23) 



The width of the distribution is proportional to l/i?o(r) 
and grows steadily as the energy scale is decreased. This 
is shown in Fig. yj for a flow near the critical point, but it 
is also true away from the critical point, and is consistent 
with the simple approximation to the RG flows discussed 
above for which the field distribution at low energy scales 
is always an exponential whose width never decreases. 

As initial conditions, we choose for convenience the 
simple exponential distribution of log-fields with the ini- 
tial Rq = 1. [Note that the strong-randomness RG equa- 
tions, Eqs.(p - ph, are invariant under a multiplication 
of all the log-couplings by any constant. Thus, although 
the initial choice Rq — 1 would appear to correspond to 
moderate randomness, we can, without loss of generality, 
use this in our study of the strong-randomness RG flows.] 
To search for a fixed point we measure l/-Ro(r) and scale 
both (3 and C by this width, defining 



/3sc — -Ro/3 



and 



Csc = -RoC- 



(24) 



(Note that Hq and thus F are not defined precisely, so we 
can not simply scale the data by F.) The scaled field dis- 
tribution is now i?sc(/3sc) — exp(— /3sc) and we can con- 
centrate on the scaled bond distribution Psc(Csc)- The 
shape of the bond distribution evolves continuously and 
its characterization is much less clear. Numerically we 
observe that for d > 2 at and near the quantum critical 
point the cluster RG always generates positively slop- 
ing {dP/d( > 0) bond distributions. This is in contrast 
to one-dimension, where the exact critical fixed point 
bond distribution Psc(Csc) = exp(— Csc) is the same sim- 
ple expoaential as the field distribution due to a duality 
relation.Q 

For d = 2 the bond distributions that are generated 
by the RG near the quantum critical point can be rea- 
sonably approximated, in the small Csc regime of inter- 
est, by a simple linear fit: Psc(Csc) — a + ^Csc- We thus 
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FIG. 1. RG evolution of the field distribution at a pu- 
tative quantum critical point. Initial conditions are a 
290 X 290 (iV = 84, 100) triangular lattice with couplings 
drawn independently from the field distribution R(P) = e~^ 
and the scaled bond distribution Psc(Csc) ~ 0.1 + 0.105(^sc, as 
described in the text. Lines are fits to the simple exponential 
form R — R(,e~^° with Rq depending on the energy scale. 
Note that the renormalized field distribution fits this form 
well for all A'^. Inset: width of the field distribution, l/i?o, 
vs. the number of remaining spins N\ the RG evolution is in 
the direction of decreasing N . The increasing width indicates 
the RG flow towards infinite randomness. The line here is a 
power-law fit that gives our estimate of the exponent tp. Note 
that this fit works well only after the rather strong initial 
transient. 




FIG. 2. RG evolution of the scaled bond distribution 
-Psc(Csc) at the apparent critical point corresponding to the 
initial conditions of Fig. hi. Note that although the intercept 
is fairly stable at Psc(O) — 0.15 in the later stages of the RG 
(this is our criterion for locating the critical point — see also 
Fig. 3), the shape of the distribution is not as stable (see also 
Fig. i ) , presumably due to transient effects from our uncorre- 
lated initial conditions. Because we do not keep bonds weaker 
than Jmin (see text), the range of the scaled bond distribution 
that we measure and plot, < i^ac < Cm = Ro\og{Q,/ J-mi-a), 
decreases as we run our RG. 



choose for initial conditions a Psc(Csc) of this form. Our 
initial lattice for all the data presented, here is a trian- 
gular lattice with periodic boundary conditions (we also 
tried others, such as square, to confirm that the results 
did not depend strongly on this arbitrary choice). Since 
we expect the stronger bonds to be shorter-ranged, we 
select the nearest-neighbor bonds (there are three such 
bonds per site) to constitute the strongest-bond part 
< Csc < Cc of the bond distribution, with Cc chosen 
so that /„ (a -I- hC,)dQ = 3, i.e. there are precisely three 
bonds per site with Q^c between zero and ^c- Then the 
next batch of the distribution, C,c < Csc < Cm, are assigned 
at random to all the second- and third-neighbor bonds 
(six more bonds per site), with Cm chosen appropriately. 
(This Cm sets our Jmin, as discussed above.) Thus our 
initial condition has nine bonds per site, corresponding 
to a coordination number of eighteen. Under renormal- 
ization, the lattice is quickly randomized, so it no longer 
resembles the initial triangular lattice, and the number 
of bonds kept per remaining cluster initially increases. 
However, since we do not keep bonds with J < Jmin, 
the number of bonds kept per cluster decreases in the 
later stages of the decimation, as fi decreases towards 
Jmin- Note that our choice of the part of the distribution 
with Csc > Cc (the tail) is only a matter of convenience; 
its details (and even its very presence) are not impor- 
tant: specifically, we have checked that by the time the 
RG reaches this energy range, most of the original bonds 
from the tail are gone, having been dominated by the 
stronger interactions that arise from the original strong 
first-neighbor bonds only. 

We first searched for a fixed point of the RG by start- 
ing from such a linear distribution of initially uncorre- 
lated bonds and monitoring the flow of the two parame- 
ters obtained by fitting the renormalized and scaled bond 
distribution to such a linear form — the intercept Psc(O) 
and the slope dPsc/dC,sc- By choosing an initial bond dis- 
tribution close to the fixed point distribution we tried to 
minimize the transients that occur as the full fixed point 
joint distribution of the fields and bonds is generated by 
the RG. However, the transients remained too strong for 
us to accurately locate a fixed point of the RG fiow in 
the plane of these two parameters: we could not fully 
stabilize this scaled bond distribution. Fig. shows an 
example of the evolution of the scaled bond distribution 
for initial conditions near what we estimate to be the crit- 
ical fixed point. The intercept stabilizes at Psc(O) = 0.15 
(see Fig. y), but the slope is much less stable, although 
it may be approaching a limit, as shown in Fig. H. 

Since we could not obtain a clear fixed point in the 
two-parameter space defined by the simple linear fit to 
the bond distribution, we instead chose as our candidates 
for critical points those that produce a scaled bond distri- 
bution whose intercept Psc(O) appeared to be stabilizing 
to a fixed point value under the action of the RG. In 
running the RG, at each decimation step the maximum 
energy term that is "integrated out" is either a bond or 
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FIG. 3. RG evolution of the intercept Psc(O) of the scaled 
bond distribution for different initial conditions: all are ini- 
tially 290 X 290 triangular lattices with independent cou- 
plings and initial field distribution R{f3) = e~ , the different 
curves correspond to various initial scaled bond distributions: 
Psc(Csc) = a + &Csc with the parameters a and b as indicated 
(see text). If the intercept saturates to a finite nonzero value 
as A'^ is decreased, this indicates that the system is critical. 
Our best estimate of the critical point has a = 0.10, 6 = 0.105 
(filled squares); this is what is used in all the other figures. 
Other parameters which we view as possibly critical are also 
indicated by filled symbols. The error estimates on the var- 
ious critical exponents include the results from all of these 
potentially critical systems. 
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FIG. 4. Evolution of the slope dPsc/d,(sc of the scaled bond 
distribution for the same set of different initial conditions as in 
Fig. H. At each energy scale, the slope is calculated by fitting a 
linear function to the corresponding distribution (as in Fig. Q) 
in the full available region < Csc < Cm- The observed shape 
of the bond distribution is only approximately linear, and 
our data for the bond distributions becomes very limited and 
noisy for small A*'. The strongly transient behavior seen here 
is presumably due to both actual transients in the shapes of 
the distributions and the reduction with decreasing A^ of the 
range, to Cm, over which the linear fit is made. 
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FIG. 5. Width of the field distribution 1/Ro vs the 
log-energy scale F for the same set of different initial con- 
ditions as in Fig. H. There is a clear initial transient in all the 
data for 1/Ro < 1.5 (see also plot of 1/Ro vs A'^ for our can- 
didate critical point in Fig. hi). After the initial transient the 
data are consistent with the expected linear behavior at crit- 
icality and with the saturation of the width in the disordered 
phase. 



a field. The ratio of the frequency of these occurrences 
is simply the intercept Psc(O) (since we have normalized 
so that the intercept of the scaled field distribution is 
Rsc{0) = 1), and is easily estimated without fitting any 
distributions by counting the number of occurrences of 
the two types of decimations as the RG runs. Stability 
of the ratio of the frequencies of the two types of deci- 
mations is a direct indication of a "balance" between the 
field part of the hamiltonian and the bond part. A de- 
tailed study, extended down to a factor of over 1000 in 
A'', shows three different types of behavior of the inter- 
cept, which we interpret as follows: If Psc(O) is clearly 
decreasing towards zero then the system is in the dis- 
ordered phase. If Psc(O) is steadily increasing then the 
system is in the ordered phase. Finally, if Psc(O) appears 
to be saturating at some value then we have a candidate 
for the critical point. Fig. 13 illustrates this. Starting 
from different initial two-dimensional lattices and differ- 
ent initial distributions, we always find that the appar- 
ently "critical" (stable) value of the intercept is in the 
range 0.1 — 0.2, using what we think are conservatively 
large uncertainties on when and where the intercept sta- 
bilizes. This range of candidate critical points, shown by 
the filled symbols in Fig. bI are used in all of our error 
estimates. 

In the simple approximation to the RG that neglects 
correlations involving fields and has Pq = -P(O) ~ l/F, 
the intercept is Psc(O) « I/-S0; in general at the criti- 
cal point Pq Ki Zq/T and the intercept is thus Psc(O) ~ 
Po/Ro = Zq/Bq. 

Fig. |] and Fig. g show the evolution under the RG of 



the field distribution and the scaled bond distribution of 
one candidate for the critical point. Since the bond dis- 
tribution is not fully stable, our scaling analysis of the 
critical flow, which we discuss next, is not as certain as 
our conclusion on the nature of the critical fixed point: 
i.e. that it is at infinite randomness. As we already men- 
tioned, direct scaling with T requires estimating the ad- 
ditional parameter Oq- To estimate the "tunneling scal- 
ing" exponent ip we therefore consider the evolution of 
the width l/i?o of the field distribution; this is shown in 
the inset of Fig. 0. It is expected that I/Rq ~ Ar-V-Zd at 
the critical fixed point since at asymptotically low energy 
scales, for which the "bare" scale Hq is not important, 
l/i?o should be proportional to F. (Fig. || shows that 
this is indeed true in the later stages of the renormal- 
ization after the initial transient.) It can be seen from 
the inset in Fig. y that during the initial transient the 
width of the field distribution grows more slowly than 
later in the renormalization. This accelerating growth of 
the width occurs for all our candidate critical points for 
d = 2, and emphasizes that the RG flow is certainly to- 
wards inflnite randomness. Fitting the later stages of the 
RG for all the initial conditions that appear consistent 
with being critical - illustrated by the filled symbols in 
Fig. H - gives exponent estimates in the rangeE3 



t/j = 0.42 ± 0.06 



(25) 



For all candidate critical points the estimated exponent tp 
is noticeably larger than the tp = d/{l + Bq) = 0.2 — 0.3 
from the simple approximation discussed above if one 
uses the 1/Bq = Psc(O) = 0.1 — 0.2 obtained from the 
apparent intercept; this indicates that correlations be- 
tween the fields and the bonds must be substantial at 
the critical fixed point. Indeed, after renormalization, 
correlations between a field and its adjacent bonds are 
easily detected: the strengths are anticorrelated, so that 
for example a cluster with a weak renormalized field is 
more likely to have strong renormalized bonds connected 
to it. 

The fractal dimension df — (pip of the critical clusters 
can be obtained directly from the RG flows at the crit- 
ical flxed point. Fig. g shows the scaling with N of the 
average magnetic moment (proportional to the number 
of bare spins) of surviving clusters. Direct flts to such 
plots for our candidate critical points give 



l.OibO.l 



(26) 



in contrast to the prediction from the simple approxima- 
tion of d/ = 0.7-0.9 . Note, however that if Sq = 4 
is obtained from our ip = 0.4 by using the simple ap- 
proximation (but ignoring the estimates of Psc(O)), the 
predicted d/ ^ 1.0 is close to the value obtained from 
the full RG. This suggests that the correlations between 
the moment and the fleld on a cluster are reasonably well 
captured by the simple approximation. 

From the scaling relation Eq.(O), the average critical 
correlations decay with the exponent 
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FIG. 6. Scaling of the average magnetic moment per clus- 
ter, Jl, with the number of remaining clusters A*'. Under the 
RG, the magnetic moment of a cluster, /i^, and local log-field 
magnitude, /3i, become strongly positively correlated, and we 
expect significant transients if our initial conditions do not 
have these correlations. The initial transient is clearly seen if 
we start with fii = 1 for all sites, or with fii = pti — an at- 
tempt to imitate the positive correlation — , but the transient 
dies off quickly as we run the RG. This transient is slightly 
suppressed if we generate initial (/3i, pi) from the joint distri- 
bution function B{pl, fj,) which is the fixed point joint distribu- 
tion in the uncorrelated-field approximation (discussed in the 
text) with B(0) = 7. We consistently find ^^ ~ Tv^o.soio.os 
implying df = 1.0 ± 0.1 and 77 = 2.0 ± 0.2. 



?7 = 2.0±0.2 



(27) 



More direct fits for the exponent (p alone can be obtained 
from plots (not shown) of the average magnetic moment 
vs. Ro, giving 



2.5 ±0.4 



(28) 



Xhe recent quantum Monte Carlo study by Pich et 
a/.til of the two-dimensional random Ising ferromagnet 
has found evidence that the width of the distribution 
of the logarithms of characteristic energies grows with 
sample size at the quantum critical point, as for one- 
dimension. They estimate ip = 0.4, in good agreement 
with what we find from the numerical RG. They also 
measured the spatial correlations G{r) —< (Tq cr^ > at 
criticality and found that the median (and hence typi- 
cal) correlation Gtyp(r) falls off faster than a power of r, 
better fit by — InGtyp ~ r^^ with tpc — 1/3, not incon- 
sistent with the scaling prediction ipc — ip- In contrast, 
the average critical correlations exhibit a power-law de- 
cay with 77 = 2, which implies that the fractal dimension 
of the critical cluster is df = 4>tp = 1, again in good 
agreement with the exponents estimated from our RG 
study. 



V. ORDERED AND DISORDERED PHASES 

We now turn to a discussion of the ordered and dis- 
ordered phases. Here and henceforth, we will denote 
the parameter that controls the difference between the 
strengths of the typical random fields and those of the 
typical random bonds in the original Hamiltonian by 
6, chosen so that the zero-temperature quantum critical 
point corresponds to S — 0, the zero-temperature disor- 
dered phase to S > and the zero-temperature ordered 
phase to S < 0. 

In our numerical RG studies, we do not know the fixed 
point accurately enough and do not have sufficient con- 
trol over initial transients to study the off-critical flows 
directly. Nevertheless, we can still obtain some infor- 
mation about the near-critical properties indirectly from 
the critical flows (just as, in conventional systems, the 
correlation length exponent, i>, is related to the decay of 
energy density correlations at the critical point). 

The effective field of a cluster is generally a product of 
some number / of the original fields divided by a product 
of (/ — 1) original interactions (both the original fields 
and original interactions need not be distinct). At the 
critical point we expect 



with the new exponent p satisfying 
p> max(<^, -). 



(29) 



(30) 



The first inequality obtains because any spin which is 
active in a cluster contributes (at least once) its origi- 
nal field to the effective field of the cluster; the second 
inequality follows from the observation that the (deci- 
mated) bonds which hold the cluster together also con- 
tribute to the effective field and must reach across the 
diameter of the cluster which is of order r~. 

An instructive way to understand the effect of devia- 
tions from criticality, is to move away from it by simply 
multiplying all the original fields by an amount 1 + S. 
Perturbatively, this would change the log-fields at scale 
r by of order fd ~ F''^. If we neglect the effects of 
d on changing the order of decimations, we obtain the 
crossover scale away from criticality 



Ta 



I'^r 



(31) 



as the scale at which the log-field changes become com- 
parable to a typical log-field or log-interaction (~ F) and 
thereby substantially alter the distribution. This can 
be justified as follows: Before the scale F^ is reached, 
S only changes the global sequence of the decimations 
but not significantly the local sequence; thus the esti- 
mate of the crossover scale is consistent. Note also that 
"chaotic" behavior under the RG flows, as occurs at, (and 
below) the critical point in classical spin glasses,cZI can- 
not occur here. This is because the classical Ising model 
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FIG. 7. Scaling of the average cluster "history" / — the 
number of the original fields whose product enters the effec- 
tive field — with the number of remaining clusters N. Under 
the RG, / and the log-field magnitude /3 become positively 
correlated, but the transient generating this correlation dies 
off quickly. By starting either from fi — 1 and the corre- 
sponding quantity for the bonds fij = 1, or from fi — Pi and 
hj = Oj, we consistently find / ~ iV-0-67±o.o7 _ ^i.34±o.i4^ 



that is equivalent to the random quantum Ising model — 
the higher dimensional generalization of the McCoy- Wu 
modelu — is purely ferromagnetic; this implies that in- 
creasing an original field on one spin in the quantum 
model must increase the last-to-be-decimated field in any 
fixed region around this spin. 

The scaling argument above yields a correlation length 

e - Cav ~ |<5r'^ with 



1 



(32) 



(p-l)V^- 

Numerically, we obtain from fitting to f ^ JSf-p^f^^ 

pV^ 1-34 ±0.14 (33) 

(Fig. 0) and hence v. We can also fit more directly for v 
using fRo ~ N~^^^'^'^\ obtaining a similar estimate 

i/= 1.07 ±0.15, (34) 

which is consistent with the boundsll3 v > 2/d = 1 and 

1 



ly < 



max(l, (j)^) — ip 



= 1.7. 



(35) 



A. Disordered Phase: Correlations and 
GrifHths-McCoy Singularities 

In the disordered phase, i5 > 0, the average spin corre- 
lations will be dominated by rare large clusters and decay 
exponentially: 



G(r) - e-^-l^ 



(36) 



The typical correlations, on the other hand, will decay 
more rapidly, as 



with 



6 



typ 



InGtyp(r) -r/^typ 



^-(i-V)'' ^ ^i-V" <• t 



(37) 



(38) 



This can be seen as follows: At the crossover scale away 
from criticality, F^, which corresponds to length scale ^, 
the ratios between typical remaining bonds and fields are 
of magnitude ln(J/ft,) ^ — F^. Two spins of interest sep- 
arated by distance r much longer than ^ will each have 
the maximum — albeit small — component of their spin on 
a remaining cluster near to them. As these clusters and 
those between them are decimated until eventually the 
two spins are contained in the same (mutual) cluster, a 
multiplicative factor of order the typical J jh ratio at the 
crossover scale will reduce the components of the spins 
on the remaining clusters — and thus on the eventual mu- 
tual cluster — for each one of the clusters at the crossover 
scale that is decimated; i.e., for each element of length of 
order ^. The result Eq.(38) follows. Note that the typical 
correlations at and near criticality have the scaling form 



InG 



typ 



r^F,y^{rli). 



(39) 



The behavior of the typical correlations is related to 
the behavior of the distributions of fields and bonds. 
The field distribution in the disordered phase has finite 
width in the limit of low energy (as in one dimension), 
i?(/3;F|(5) -^ i?o((5)e--«o(5)/3 for F > F^, while the bonds 
continue to become weaker. The limiting width is of or- 
der 

\IRa{S) ^ F, ^ S-^'' (40) 

which diverges as (5 -^ O.ErtJ As in one-dimension, this 
gives rise to a disordered "Griffiths" phase with contin- 
uously variable dynamical exponent z that relates the 
scales of energy and length — the typical fields and spac- 
ing of surviving clusters — via 



9. 



(41) 



with z diverging as 



z - b-'^'' (42) 

for (5 ^ 0, consistent with the numerical data of Ref.llil. 
In the disordered phase the distribution of log-fields 
tends to a simple exponential ffujpji at low energies with 
the limiting width l/i?o ~ z/d.oQ Concomitantly, there 
is a continuously variable power law singularity in the av- 
erage ground-state magnetization per spin in an ordering 
field H: 



with the exponent of the logarithmic factor not deter- 
mined from these simple arguments. This gives the lead- 
ing low-i/ behavior near the critical point where z > d. 
For larger S where z < d, this instead gives a singular 
correction to an analytic M{H). The low temperature 
zero-field susceptibility likewise diverges for z > d as \/T 
raised to a continuously varying power that is less than 
unity — a weaker than Curie law divergence. 

For small (positive) H and small |(5|, the magnetization 
has a scaling form: 

M[H,5) - |lni/|'^-''/'^S(C(5|lni/|i/'^'') (44) 

with both C and the cutoff scale implicit in the h\H 
representing non-universal corrections to scaling. When 
u is large and positive, the scaling function S(u) ^ 
exp(— u'^'') times a power of u. 

The clusters surviving to low energies in the disordered 
phase are rare large strongly coupled ferromagnetic clus- 
ters which exist even in the disordered phase. The most 
probable way for such a region to occur is (it appears) 
for there to be, at the crossover scale F^, a connected set 
of n clusters each with linear size ~ ^ with somewhat 
anomalously weak fields on them and somewhat anoma- 
lously strong effective couplings between them; these will 
make them join together at lower energies into the rare 
large cluster of interest. Note that because the system 
has strong randomness, there is no requirement that this 
large cluster has a compact geometry. Indeed, for large n 
we expect that the most probable such clusters will have, 
on scales larger than ^, a geometry similar to incipient 
classical percolation clusters, as the requirements for a 
large rare cluster to form out of the scale ^ sub-clusters 
is qualitatively like that for a large cluster to form out of 
the small scale objects in conventional percolation. [See 
more detailed discussion in next sub-section.] Since the 
probabilities of the occurrence at scale ^ of each of the n 
such clusters and each of the (~ n) corresponding cou- 
plings are determined primarily by local properties and 
are hence roughly independent, the probability of the 
rare large cluster is of order l/a" with a > 1; a is in- 
dependent of 5 for small 5, since the needed clusters and 
bonds are not very atypical at scale ^. The field on such 
a cluster will be of order the energy scale for crossover 
away from criticality, multiplied by of order n factors, 
each less than unity with logarithm of order F^: 



r!o(e-'^^^) 



(45) 



with a independent of 5 near the critical point. The 
probability of a cluster with an anomalously small field 
h is thus approximately 



p{h) 



^-{\na)\n{Q.a/h)/ars 



(46) 



M r-. H'^/''\\nH\ 



(43) 



If one finds the least rare such clusters, i.e. with (lnQ;)/a 
as small as possible, then Ro{d) — {\na)/aTs and these 
rare clusters give rise to the singular scaling in the Grif- 
fiths phase quoted above. As stated earlier the dynamic 
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exponent z is given by the asymptotic low energy value 
of d/RQ. [Note that we have ignored here differences 
between In fig and In ^g = In fig ^ r* since these will 
not affect the dominant behavior in the limit of interest. 
Also note that since \hih\ ~ n, a factor p{h) controls 
the asymptotic behavior of the distribution of the | \nh\- 
variable.] 

The RG generates the low-energy tail r-^ q-^o{S)I3 qJ 
the log-field distribution on scales L ^ ^ precisely be- 
cause the rare large clusters discussed above are not too 
rare. Note that an exponential tail of the distribution of 
the cluster effective fields is generated by the RG almost 
immediately and some rare arbitrarily large clusters are 
present at any energy scale, as long as there is no infinite 
cluster that would consume them. In particular, this tail 
is also generated at the critical point, but it continues to 
become longer as the energy scale decreases. 



B. Ordered Phase: Percolation and Finite 
Temperature Ordering 

The behavior in the ordered phase differs strikingly 
from that in one dimension. At a finite energy scale 



i^oo — f^oe 



-r„ 



noeM-K\5\~''"'] 



(47) 



a single infinite cluster (with zero transverse field) devel- 
ops. Finite clusters, some of which will join the infinite 
cluster at lower energies, coexist with it. The sponta- 
neous magnetization is proportional to the number of 
sites in the infinite cluster at fi ^ yielding 



with 



Mo ^ \S\^ 



(3 = vr]/2 = v{d--(l)'ip). 



(48) 



(49) 



This implies that the scaling function in Eq.(^^ has the 
asymptotic form in the ordered phase S(u —^ — cxj) ^ 

As a consequence of the infinite cluster development, in 
the ferromagnetically ordered phase there is, in contrast 
to one-dimension, an interfacial energy density — albeit 
exponentially small — for a domain wall for any S < 0. 
This 'surface tension' scales as 



e-'n. 



(50) 



— simply the energy needed to cut the infinite cluster in 
half. The transition temperature for destruction of the 
long range order by thermal fluctuations will similarly 
be determined by the strength of the bonds that hold 
together the infinite cluster; thus 






(51) 



We note that the arguments given below imply that the 
width of the classical finite temperature critical region as 



the temperature is reduced at fixed negative S, will be of 
order T^. 

The behavior near the energy scale r2oo((5) at which 
the infinite cluster forms has aspects similar to con- 
ventional classical percolation (in contrast to the zero- 
temperature quantum percolation transition at S = 0). 
At the crossover scale Ts = In^flo/ils) away from crit- 
icality, the bonds will start to dominate over the ran- 
dom fields. At this scale the sizes of the clusters and 
the lengths of the bonds will typically be of order the 



crossover length scale £, 



-i/V-. 



indeed, bonds much 



longer than this will be exponentially rare. Between the 
log-energy scales Ts and Too ~ ln(rio/^oo), at which the 
infinite cluster forms, most of the decimations will be of 
bonds, resulting in the joining together of clusters. The 
process of decimation of bonds (and occasional clusters) 
will continue with larger and larger clusters forming until 
the percolation scale T^o- Note that T^o will be a fixed 
(order-one) multiple of the somewhat arbitrarily defined 
crossover scale Tg , but the corresponding physical energy 
fioo is actually exponentially smaller than the energy ^s 
(for large Ts). 

Unlike the case at the quantum critical point, the pro- 
cess by which the large clusters are joined together as the 
energy scale is decreased near the percolation scale is ba- 
sically local. The key feature of this locality, which occurs 
when r is in the range F^ < F < Foo, is that when F is 
changed by a small amount, whether, say, a large cluster 
A will become joined to a large cluster B, and whether 
the same cluster A will become joined to another large 
cluster C, are roughly independent events, each only de- 
pending on the smaller clusters and bonds - which have 
typical length scale ^ and log-energy scale F^ - in the 
vicinity of the respective potential connections. As the 
percolation scale is approached, we expect that this inde- 
pendence will become more and more pronounced as the 
important connections that make the large clusters grow 
become further and further apart. On the basis of this 
argument, we conjecture that the percolation that occurs 
in the cluster RG at scale Fqo is in the universality class 
of classical percolation, with (Too — T) playing the role of 
(Pc — p) in classical percolation. [Note that the width of 
the log-bond distribution at scale F^ is of order F^ and 
will remain so at all scales F^ < F < Too]. 

The nature of the percolation process at the scale Too 
controls the critical behavior associated with the finite- 
temperature ordering transition at T^ '^ fioo over a sub- 
stantial region of the T—S plane. On the logarithmic tem- 
perature scale. Ft = ln(f2o/T), the finite-temperature 
spin-spin correlation length, S.T, is simply the character- 
istic length scale of the distribution of cluster diameters 
at scale Fy: clusters that exist at this scale will each 
consist of sets of well-correlated (active) spins, while the 
correlations between these clusters will be destroyed by 
the thermal fluctuations. On temperature scales above 
i^s {Tt < Ts), £,T will be given by the quantum criti- 
cal length scale at log-energy scale Ft. But on scales 
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between fls and ^oc {Ts < Ft < Too), the percolation 
process will cause the correlation length to diverge as 



CT'-e 



r^ 



r. 



Tn 



(52) 



where Vp is the classical percolation correlation length 
exponent — with Vp = 4/3 in two dimensions — and 
^ ^ [—5)~'^ is the correlation length associated with the 
quantum critical point (i.e., the characteristic length at 
the crossover scale r^-). 

As the critical temperature is approached, the RG ap- 
proximation will eventually break down at any non-zero 
d. The clean separation of bonds into 'strong' for those 
stronger than Ft, and 'weak' for those weaker than F^ 
will not hold for the 'marginal' bonds whose strength is 
of order T, which correspond to those whose log-strength 
is Ft ± 0(1). This implies an 0(1) multiplicative uncer- 
tainty in the proportionality between Tc and Oqo and it 
also implies that the percolation-dominated form of the 
critical behavior of Eq. (p2) breaks down when Fqo — F^ 
is of order one. Closer to the finite-temperature criti- 
cal point, the behavior will be dominated by the thermal 
fluctuations of the marginal bonds that link very large 
almost percolating clusters. This will make the critical 
behavior cross over to that of the conventional classical 
d-dimensional random bond Ising universality class with 
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where vj is the classical random Ising correlation length 
exponent, equal to one in two dimensions (assuming the 
interactions are not frustrated). 

As a function of temperature, this double crossover 
in the critical behavior will be particularly hard to ob- 
serve due to the logarithmic temperature scale which 
makes the crossover energy scale VLg exponentially larger 
than Tc for small (negative) 5. But if the tempera- 
ture is held fixed and extremely small — i.e. Ft ^ 1 — 
then the crossovers can be seen more readily by de- 
creasing the relative strengths of the random fields that 
are parametrized by 5. The critical value, 5c, is of or- 
der — F^ . As the random fields are reduced from 

5 of order one until 6 ~ +Tj, ' ~ |i5c|, the ther- 
mal effects will be negligible and ^t will diverge with 
the zero-temperature quantum, exponent that we have 
denoted simply v. As 5 is further decreased through 

zero until i5 — i5c ~ F^ ' "^ ~ |i5c|/Ft, the interme- 
diate classical-percolation dominated critical behavior as 
in Eq. ( |52| ) will obtain with the exponent Vp. Finally, for 

5 — 5c^ F^ ~ ' the classical random Ising critical ex- 
ponent vi will control the divergence of ^t. Note that in 
the limit of asymptotically small l/Fy, all three of these 
regimes will become very broad on a ln(((5c — 5) /5c) plot. 



C. Ordered Phase: Singularities 

We now turn to properties of the ordered phase on 
energy and temperature scales much lower than the or- 
dering temperature Tc. The low-energy properties of the 
ordered phase for d > 1 will not have the strong power- 
law Griffiths singularities found in the disordered phase 
and the one-dimensional ordered phase. This can be seen 
in the RG language which naturally incorporates the role 
of rare anomalous regions. If we continue the RG much 
below VLoc, i.e., beyond the formation of the infinite clus- 
ter, we will find that the remaining finite clusters are 
almost always connected only to the macroscopic clus- 
ter and at lower energy scales almost all of these will 
either be decimated or will join the infinite cluster; fewer 
and fewer will join together to make larger finite clus- 
ters. Since the finite clusters and the bonds connecting 
them to the macro-cluster are decimated independently, 
no new clusters or bonds will be generated, and the low- 
energy tails of the distributions remain essentially the 
same as they were just below Q,rx, (when the widths of 
both the log-field and the log-bond distributions are of 
order ~ \n{no/n^) ~ F5 ~ \5\~^''). 

As in the disordered phase, the pre-formed tails of the 
field and bond distributions represent rare large regions 
responsible for the low-energy excitations in the system. 
The role of these rare fluctuations, however, is very dif- 
ferent in the ordered phase; although they still make the 
system gapless, they do not dominate all the low-energy 
properties (such as, for example, the response to a small 
ordering magnetic field). In contrast to the disordered 
phase, the dominant rare regions in the ordered phase 
are indeed very rare and do not produce a power-law sin- 
gularity in the density of states at zero energy. This can 
be seen by analyzing the probability that a cluster with a 
very small effective field h survives down to energy scale 
Vt ^ h <^ r^oo. We can consider such a surviving ferro- 
magnetic cluster to be composed of n subclusters (each 
of diameter ~ £) with h ~ J7oe~™'"'', as in the disordered 
phase. But this cluster must be isolated very effectively 
from the rest of the system — with effective coupling link- 
ing it to the infinite cluster of order h or weaker. In the 
disordered phase, the typical length of a bond with effec- 
tive coupling I In J| ~ nVg is L ~ n^. Thus to achieve suf- 
ficient isolation, the disordered region around the droplet 
must have a linear size of order r^ n^. The probability of 
such a rare region is very small — of order 1/a" — so that 
the generic low-energy tail that the RG can generate is 
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(54) 



very different from the ^ p-c\inh\ ^^^^ in the disordered 
phase and the similar tail in the one-dimensional ordered 
phase. It is also strikingly different from an even longer 
tail, ~ g-c|in/t| ^ -j-j-^g^-j- Qccurs in the cpdered phase of 

the dilute quantum Ising system of Ref.ll3 which differs 
from ours by some fraction of the initial Jij 's being zero, 
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so perfectly isolated clusters can form even in the ordered 
phase. The origin of the difference between these cases is 
easy to understand: in one-dimension, the length and the 
volume of an isolating region are the same, while in the 
dilute case, with a delta-function weight at zero coupling, 
a droplet of size n can be completely isolated by a "dis- 
ordered" region of volume n^'^~^'''^ — just a surrounding 
surface of missing bonds. 



the quantum critical point is approached, these rare clus- 
ters and their couplings will effectively act like a strongly 
random system which we expect will dominate the behav- 
ior and cause the whole system to be driven to strong (but 
not necessarily infinite) randomness sufficiently close to 
the critical point - however weak the original random- 
ness. This intriguing possibility clearly merits further 
investigation. 



VI. HIGHER DIMENSIONS AND RELATED 

SYSTEMS 



B. Other Quantum Transitions writh Discrete Broken 
Symmetries 



A. Higher Dimensional Random Ferromagnetic 
Quantum Ising Models 

So far, we have presented detailed results only for two 
dimensions, although the general scaling picture, expo- 
nent equalities, behavior of correlation functions, etc., 
should be qualitatively the same for strong randomness 
in any dimension (d > 1) for which there is a stable 
infinite-randomness critical fixed point. Our numerical 
studies in three dimensions are sufficient to indicate that 
the infinite-randomness fixed point is stable, although 
they are not thorough enough to yield reliable estimates 
of exponents and their uncertainties. Since in both two 
and three dimensions, weak randomness is a relevant per- 
turbation away from the pure fixed point, we expect the 
same strong-randomness-dominated critical behavior to 
occur even for arbitrarily weak randomness. 

The situation in higher dimensions - d > 4 - is far 
more uncertain. It is not clear at this point whether or 
not the direction of the RG flow at strong randomness 
reverses for d sufficiently large. 

For weak disorder, it would appear that the situation is 
clearer: the Harris criterion would seem to indicate that 
for d > A weak randomness is irrelevant. But one must be 
very careful. There are other situations known in which 
weak randomness formally appears to be irrelevant but 
for which exponentially rare regions change the behavior 
for arbitrarily weak randomness.ta We strongly suspect, 
as argued below, that this will be the case here. In gen- 
eral, Griffiths singularities and other strong-randomness- 
like effects will start to appear when the random quantum 
Ising system is close enough to the quantum transition 
that the distribution of J's and the distribution of ft,'s 
overlap in the sense that for some values of {h, J) in the 
support of these distributions, a pure system would be in 
the ordered phase, while for other values of {h, J), a pure 
system would be in the disordered phase. This implies 
that arbitrarily large rare regions will exist that act as if 
they were in the opposite phase than the full system is. In 
particular, in the disordered phase sufficiently close to the 
critical point, strongly correlated clusters will exist with 
broadly distributed effective fields and effective interac- 
tions between them which are broadly distributed and 
typically decay exponentially with their separation. As 



As mentioned in the Introduction, the infinite random- 
ness critical fixed points found here control more than 
just Ising ferromagnetic quantum transitions. In par- 
ticular, as pointed outrfor the one-dimensional case by 
Senthil and Majumdaifl, Potts models or any random 
quantum systems with continuous (second order) transi- 
tions at which a discrete symmetry of a non-conserved 
order parameter is broken, will have the same critical be- 
havior as the Ising case, with the extra degrees of freedom 
just "going along for the ride" on the basic geometrical 
transition. 

This holds even for systems which are frustrated on 
small length scales, such as quantum Ising spin glasses. 
Because of flow to the infinite randomness fixed point, 
the frustration will become irrelevant at low energies at 
the critical point, since in any frustrated loop the weak- 
est interaction will be infinitely weaker than the others, 
so can be ignored. The primary changes here concern the 
coupling to a uniform magnetic field in the z-direction, 
and the behavior at nonzero temperature in the ordered 
phase. Because the uniform field is not an ordering field 
for the spin glass, the magnetic moment of the large clus- 
ters will be random in sign, scaling as the square root of 
the number of active spins on the cluster. At the critical 
point this will change the M{H) scaling, yielding 



MsG-- |lni/|5<^-'^/^ 



(55) 



in contrast to Eq.(|21|) for the ferromagnetic case. In the 
disordered phase, M will scale as the same power of H 
as in the ferromagnetic case with only the logarithmic 
prefactors modified. In the ordered phase, the behavior of 
the nonzero-temperature, long-distance correlations will 
cross over to classical spin glass behavior at and below 
a temperature of order f2(x>. For d = 2, true long-range 
spin-glass order will be present only at zero temperature, 
because the lower critical dimension for the classical spin 
glass is always more than two. 



C. Random Quantum XY And Heisenberg 
Antiferromagnets 

The simplest example of an infinite randomness quan- 
tum fixed point occurs for one-dimensional random 
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Heisenberg (or XY) spin chains. In the corresponding 
phase, the "random singlet phase" , each spin is paired in 
a singlet with one other spin, usually one close by, but 
a small fraction of the spins are paired very weakly with 
spins far away. The RG analysis oLthis system, first car- 
ried out by Ma, Dasgupta and HuEl and then more fully 
by one of ug3, is a simpler version of that used in the 
present paper. 

A similar RG analysis was carried out for two- and 
three-diaiensional random antiferromagnets by Bhatt 
and Leecl over a substantial range of energy scales, in 
particular including those relevant for experiments on the 
insulating phase of phosphorus-doped silicon. This inves- 
tigation has been extended by two of ust3 to the strong- 
randomness limit. We have found that for d >2, in con- 
trast to one-dimension, the infinite-randomness random- 
singlet fixed point of random Heisenberg or XY quantum 
antiferromagnets is unstable towards a state with finite 
randomness and, presumably, more conventional scaling; 
this state includes both antiferromagnetic and ferromag- 
netic effective interactions, and involves clusters with mo- 
ments much larger than those of the single spins that 
dominate the low energy behavior in the one-dimensional 
case. 



in two dimensions, first order transitions are not pos- 
sible in random systems, this class should include all 
discrete-symmetry-breaking transitions providing there 
are no conservation laws which alter the quantum dynam- 
ics in an essential way (for an example of an Ising case 
with a conserved order parameter, seeo). The nature of 
the discrete symmetry breaking quantum transitions we 
have studied is controlled by a novel type of percolation 
— rather surprising given the intrinsic quantum nature of 
the underlying models. As the rules of this percolation 
process are asymptotically classical (although the pro- 
cess is not conventional percolation) , one might hope that 
conformal field theory approaches which take advantage 
of the two-dimensional (rather than "2 + 1" dimensional) 
structure could perhaps be used to obtain analytic re- 
sults for some of the properties of such two-dimensional 
random quantum systems. 

We thank Ravin Bhatt, Kedar Damle, Matthew Hast- 
ings and Peter Young for helpful discussions. Support for 
this work was provided by the National Science Foun- 
dation through grants DMR-9400362, -9630064 and - 
9802468 and via Harvard's MRSEC. 



VII. SUMMARY 

In summary, we have studied random quantum Ising 
ferromagnets using an energy space cluster RG which be- 
comes exact for strong randomness. Using the structure 
of the RG, we presented a scaling picture of the behav- 
ior near an infinite-randomness quantum critical fixed 
point that can occur. Near to this fixed point - corre- 
sponding to low energy scales near the zero-temperature 
quantum phase transition - the RG yields asymptotically 
exact results. We have implemented the RG numeri- 
cally, primarily in two dimensions, and found that the 
critical behavior is indeed controlled by such an infinite- 
randomness fixed point, as in one-dimension. We esti- 
mated numerically the corresponding critical exponents 
in two-dimensions, and discussed the properties of the 
disordered and ordered phases. In the disordered phase 
we found that rare anomalously strongly coupled ferro- 
magnetic clusters — in the RG language, a low-energy tail 
of the cluster field distribution generated by the decima- 
tion procedure — dominate the low-energy behavior and 
cause power-law Griffiths-McCoy singularities near the 
phase transition. In the ordered phase for d > 1, on the 
other hand, the Griffiths singularities are much weaker, 
and do not produce divergences in thermodynamic quan- 
tities; the low-energy density of states they produce van- 
ishes faster than any power of the energy. 

The universality class controlled by the infinite ran- 
domness quantum Ising critical fixed point is very broad; 
it includes all continuous quantum transitions in random 
systems at which a discrete symmetry is broken; since 
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